Initial temperature and EoS of quark matter from direct photons 
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Abstract 

The time evolution of the quark gluon plasma created in gold-gold collisions 
of the Relativistic Heavy Ion Collider (RHIC) can be described by hydro- 
dynamical models. Distribution of hadrons reflects the freeze-out state of 
the matter. To investigate the time evolution one needs to analyze pene- 
trating probes, such as direct photon spectra. Distributions of low energy 
photons was published in 2010 by PHENIX. In this paper we analyze a 3+1 
dimensional solution of relativistic hydrodynamics and calculate momentum 
distribution of direct photons. Using earlier flts of this model to hadronic 
spectra, we compare photon calculations to measurements and flnd that the 
initial temperature of the center of the flreball is at least 519±12 MeV, while 
for the equation of state we get = 0.36 ± 0.02. 



1 Perfect fluid hydrodynamics 

In the last several years it has been revealed that the matter produced in the 
collisions of the Relativistic Heavy Ion Collider (RHIC) is a nearly perfect 
fluid [1], i.e. it can be described with perfect fluid hydrodynamics. There 
was a long search for exact hydrodynamic models (solutions of the partial 
differential equations of hydrodynamics) and several models proved to be 
applicable. There are 1+3 dimensional models, as well as relativistic models 
- but no 1+3 dimensional and relativistic model has been tested yet. In this 
paper we extract observables from the relativistic, ellipsoidally symmetric 
model of ref. j2]. Hadronic observables were calculated in ref. [3]. Here we 
calculate momentum distribution and elliptic flow of direct photons. 

Perfect fluid hydrodynamics is based on local conservation of entropy 
or number density (n), energy-momentum density (chosen as T^^^ = (e + 
p)u^u'^ — pg^'^) for perfect fluids. Here e is energy density, p is pressure and 
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g^"^ is the metric tensor, diag(— 1, 1, 1, 1). The conservation equations are 
closed by the equation of state, which gives the relationship between e and p. 
Typically e = npis chosen. Note that exact, analytic result for hydrodynamic 
models is, that the hadronic observables do not depend on the initial state 
or the equation of state separately, just through the final state jl]. Thus if 
we fix the final state from the data, and determine initial state parameters 
from direct photon spectra. 

Even though many solve the above equations numerically, there are only 
a few exact solutions for these equations. See a small review of solutions in 
ref. [3]. As far as we know, until now there was only one 1+3 dimensional rel- 
ativistic solutions investigated: the solution in ref. [2]. Hadronic observables 
from this solution were computed and compared to data in ref. [3]. Present 
paper calculates thermal photon observables from this realistic 1+3 dimen- 
sional model and compares them to data for the first time. Our method is 
different from numerical calculations: here one can determine the best values 
of the parameters of the solution by fitting the analytic model results to data. 

2 The analyzed solution 

The analyzed solution, as described in refs. |2l[3] assumes self-similarity and 
ellipsoidal symmetry. The ellipsoidal symmetry means that at a given proper 
time the thermodynamical quantities are constant on the surface of expand- 
ing ellipsoids. The ellipsoids are given by constant values of the scale vari- 
able s = rl/X{tf + rllY{tf + rl/Z{tf, where X{t) = XqI, Y{t) = Yot, 
and Z{t) = ZqI are time proportional scale parameters (axes of the s = 1 
ellipsoid), and spatial coordinates are r^., r^, and r^. The velocity-field is a 
3D Hubble-type expansion, u^{x) = x^/t. This means that the solution is 
accelerationless. Note that similarly to ref. [3], we use transverse expansion 
Ut and eccentricity e instead of x and y direction expansion rates Xq and 

The temperature is T(x) = Tq (tq/t)'^^^ iy(s)^^, the number density n{x) = 
Uq (tq/t)^ with r being the proper time and z/(s) an arbitrary function 
of s. We choose tq to be the time of the freeze-out, thus Tq is the central 
freeze-out temperature. The function z/(s) is chosen as z/(s) = exp(— 6s/2), 
where b = is the temperature gradient. If the fireball is the hottest in 

the center, then b < 0. Note that for the momentum distribution of direct 
photons, we don't need any kind of density just the temperature distribution. 

The picture used in hydro models is that the pre freeze-out (FO) medium 
is described by hydrodynamics, and the post freeze-out medium is that of 
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observed hadrons. The hadronic observables can be extracted from the so- 
lution via the phase-space distribution at the FO. This will correspond to 
the hadronic final state or source distribution S{x,p). See details about this 
topic in ref. j3]. It is important to see that the same final state can be 
achieved with different equations of state or initial conditions However, 
as we will see below, the source function of photons is sensitive to the whole 
time evolution, thus both to initial conditions and equation of state as well. 
For the source function of photon creation we have then 

S{x,p)d'x = Ar ^f^"^ d'x (1) 

exp {p^u^'{x)/T{x)} - 1 

where A/" = g/{2'n-)^ (with degeneracy factor g), p^u'^ is the energy of the 
photon in the co-moving system (from the Cooper- Frye prefactor). We will 
use a second order saddle-point approximation. In this approximation the 
point of maximal emittivity and Gaussian source widths can be calculated. 



3 Thermal photon observables 

The invariant single particle momentum distribution can then be calculated 
from the source function by integrating on space-time. Our calculated quan- 
tity will be Ni{pt), which is the invariant single particle momentum distribu- 
tion taken at midrapidity and integrated on the azimuthal angle (f, similarly 
to ref. The result on the invariant transverse momentum distribution 
depends on the initial and final times. We introduce the variable ^ = t/TQ 
(with To as the freeze-out time), then in terms of ^ the time integration goes 
from i to 1, if « = tj/ro. The result is: 
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where we introduced the auxihary quantities Px = i^/ yK — — Kh/X^y and 

Py, pz similarly, while A = 1 — [px + Py)/4: and B = [px— Py)/4:. Furthermore, 
ao, ai are the Taylor-coefficients of the ffist two modified Bessel functions: 
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As the coefficients are stronly decreasing, in real calculations we can restrict 
ourselves to use only first two of them, i.e. we can make the approximation 
of Io{x) = X and Ii{x) = x^/2. 

We can also calculate the elliptic fiow, the second Fourier coefficient of the 
invariant single particle momentum distribution, similarly to ref. [3]. Result 
is: 
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Here the I2 modified Bessel function had to be introduced as well, with 
coefficients of aa = (0,0,1/8,0,1/96,0,1/3072,0, ...). 



4 Comparison to PHENIX data 

Hadronic data were described with this model in ref. [3], determining freeze- 
out (final state) parameters (expansion rates, freeze-out proper-time and 
freeze-out temperature). We use the parameters of the hadronic fit and leave 
only the remaining as free parameters. The free parameters will be k (the 
equation of state parameter) and ti, the initial time of the evolution. 



4 



Parameter 




Value 


Type 


Central FO temperatre 


To 


204 MeV 


fixed 


FO proper-time 


To 


7.7 fm/c 


fixed 


Eccentricity 


€ 


0.34 


fixed 


Transverse expansion 


uf/b 


-0.34 


fixed 


Longitudinal expansion 


zlih 


-1.6 


fixed 


Equation of State 


K 


7.7 ±0.7 


free 


Initial time 


ti 


- 0.7 fm/c 


free 



Table 1: The first five fit parametes were taken from the hadronic fits of 
ref. 13] (see more details about the parameters therein). The EoS parameter k 
was fitted, while for the initial time we determined an interval of acceptability 
(with 95% confidence). The of the fit was 7.0 here, the number of points 
was 5, so with two fitted parameters our confidence level is 7.2% 

We use direct photon data from the PHENIX Collaboration measured 
in ^/sN^ =200 GeV Au+Au collisions. Note that we also utilized a normal- 
izing factor to describe the data. Fit parameters are summarized in table [T] 
The fit itself is shown on fig. [TJ The equation of state result is k = 7.7 ± 0.7, 
or alternatively, using n = 1/c^: 

Cs = 0.36 ± 0.02 (6) 

which is in nice agreement with both lattice QCD calculations |B] and exper- 
imental results from hadronic data [31H]. Note however that the spectrum 
is not very sensitive to the initial time as in early times the thermal pho- 
ton emission is not in the region of the experimental data. Our model does 
not contain acceleration, it is a Hubble-flow type of model, but the initial 
acceleration does not play a large role in the thermal photon spectrum be- 
cause of this insensitivity on the initial time. We determined an "interval of 
acceptability" for ti. The maximum value for ti within 95% probability is 
0.7 fm/c. This can then be used to determine a lower bound for the initial 
temperature. Thus the initial temperature of the fireball (in its center) is: 

Ti > 519 ± 12MeV (7) 

at 0.7 fm/c. The uncertainty comes from the uncertainty of k. This is in 
accordance with other models' 300 — 600 MeV result [5]. 
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Figure 1: Invariant transverse momentum of direct photons from our hydro 
model. The model validity goes until roughly 3 GeV in transverse momentum. 
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